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Abstract 

C| The adsorbate dynamics during simultaneous action of desorption and lateral 

Q-i adsorbate diffusion is studied in a simple lattice-gas model by kinetic Monte Carlo 

simulations. It is found that the action of the coverage-conserving diffusion process 
during the course of the desorption has two distinct, competing effects: a general 
acceleration of the desorption process, and a coarsening of the adsorbate configura- 
tion through Ostwald ripening. The balance between these two effects is governed 
by the structure of the adsorbate layer at the beginning of the desorption process. 
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1 Introduction 

Adsorption and desorption of submonolayer structures are nonequilibrium processes of 
significant scientific and technological importance. Examples can be observed, for in- 
stance in systems that undergo underpotential deposition pQ. Here we consider the 
influence of lateral adsorbate diffusion on the time evolution of the adsorbate config- 
uration during desorption of a submonolayer from a single-crystal surface. While the 
adsorbate coverage is obviously reduced by the desorption process, it is conserved by 
the process of lateral diffusion. Thus, whereas desorption causes shrinkage of adsorbate 
clusters of all sizes, lateral diffusion may lead to a coarsening of the adsorbate config- 
uration through Ostwald ripening [2], [3]- The interplay of these two processes makes 
the resulting adsorbate dynamics quite complex. A likely experimental realization was 
observed by He and Borguet 0] during gold-cluster dissolution on a Au(lll) electrode. 
In these experiments, gold atoms were released during a short positive potential pulse 
onto the substrate, where they quickly nucleated and formed islands. After the end of 
the pulse, small islands tended to decay quickly, while large islands initially continued 
to grow before they also eventually decayed. 

In the present paper we study such a simultaneous process of desorption and lateral 
diffusion in a simple lattice-gas model by kinetic Monte Carlo simulations. The rest 
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of the paper is organized as follows. The lattice-gas model and the kinetic Monte 
Carlo algorithm are introduced in Sec. [21 the different initial configurations used for the 
desorption process are discussed in Sec. O our numerical results on desorption without 
and with lateral diffusion are presented in Sec. [U and a brief discussion and conclusions 
are given in Sec. El 



Our aim is to study how the cluster size distribution of the initial configuration influences 
the effects of lateral diffusion on the dynamics of desorption from a submonolayer of 
adsorbed atoms or molecules. The simplest kind of model for submonolayer adsorption, 
desorption, and diffusion is a lattice-gas model with only nearest-neighbor interactions. 
The grand-canonical effective Hamiltonian defining this model is 



where the first sum runs over all nearest-neighbor pairs of lattice sites, and <ft > 
denotes an attractive nearest-neighbor interaction suitable for describing underpotential 
deposition. Here, Cj is the occupation variable at site i, with a = if empty and Cj = 1 
if occupied. The second sum runs over all sites, and [i is the electrochemical potential. 
The value (jlq = — 2(f) for a square lattice corresponds to equilibrium coexistence between 
a low-coverage and a high-coverage phase at a discontinuous phase transition. This 
Hamiltonian is equivalent to a nearest-neighbor Ising spin model [5]. 

The simulations reported here were performed on a 256 x 256 square lattice with 
periodic boundary conditions at a temperature T = 0.8T C with T c the exact critical 
temperature of the Ising lattice-gas model [6]. Energy and temperature units are chosen 
such that Boltzmann's constant ks = 1. The surface coverage (particle concentration 
or fraction of occupied sites on the lattice) is given by 9 = N p /N, where N p is the total 
number of occupied sites and N is the total number of sites on the lattice. Each time 
we apply an elementary step, the system energy changes from E A to E&. To accomplish 
this, the system has to overcome a local energy maximum En. We define the height of 
the corresponding energy barrier A by the Butler- Volmer type expression [5], 



2 Model and Methods 
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The energy difference between Eh and Ea is 

AE = E U -E A 
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We then define the transition rate as 



R a ^b = exp (- (A/T)) exp (- (E B - E A ) /2T) . (4) 

We consider two kinds of elementary steps: (1) an adsorption/desorption step, and 
(2) a diffusion step to a nearest-neighbor site. We apply a rejection-free (n-fold way) 
algorithm EJ EJ [10] , in which we update the clock after every accepted configuration 
change. The energy barrier for an adsorption/desorption step is A ads / dcs , and for dif- 
fusion it is A d . The absolute values of these energy barriers are not important for our 
simulations (they are important if we want to relate the simulation time to physical 
time). Our simulations are controlled by the ratio between adsorption/desorption and 
diffusion steps. This is achieved by fixing A ads / des at a constant value throughout all 
simulations and varying A d . Attachment between adjacent adparticles is instantaneous 
in this model. 

During the simulation, we sample the time development of the coverage 9. At 
certain coverages we calculate the cluster-size distribution using the Hoshen-Kopelman 
algorithm [llj. Using this algorithm we calculate the number density p(s) which is 
the probability of finding on a randomly chosen site the center of a cluster of size s. 
During the simulation we also measure the correlation length, which for the square- 
lattice attractive lattice-gas model can be calculated as [121 02] 

where S is the total number of broken bonds at a given coverage 9. 

We start from an empty lattice and equilibrate at negative chemical potential (/i — 
/j-o < 0) to achieve a very low coverage before switching to fj, — fio > for adsorption 
until a coverage cutoff ^cutoff is reached. The configuration at # C utoff then becomes the 
initial configuration for the desorption processes with #; n it = Cutoff • The coverage cutoff 
is chosen to be 9- m \t = 0.35, well below the random percolation threshold. We prepare 
four classes of initial configurations by applying four different chemical potentials during 
the adsorption stage: \i — /io = 0.4, 1.0, 2.0, and 10.0. Figure [T] shows typical snapshots 
of the different initial configurations. 

For each run of the simulation, we performed adsorption until 9 = # C utoff> immedi- 
ately followed by desorption until a low-coverage equilibrium was reestablished. The 
simulation time t = is defined as the beginning of the desorption stage. During the ad- 
sorption stage, we fixed the adsorption/desorption barrier at A adS/ / des = 15 and turned 
off the diffusion by setting A d = 150. 

As mentioned above, we varied the chemical potential for the adsorption stage to 
obtain different initial configurations for the desorption process. For the desorption 
stage we fixed the chemical potential at /i — /ip = —0.4 and the adsorption/desorption 
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barrier at A adS/ / des = 15. To study the effects of diffusion, we performed two different 
sets of runs for the desorption stage; one without diffusion (A d = 150) and the other 
with relatively fast diffusion (A<j = 8). During each desorption run we measured 0(t), 
£(0), and p(s). We averaged the results of 100 independent runs for each simulation 
setup. 

The results of the desorption simulations are displayed in the figures, showing to- 
gether the results without and with diffusion. We have four sets of results, corresponding 
to the four different classes of initial configurations resulting from the different chemical 
potentials used during the adsorption stage. By comparing the results for the simula- 
tions without and with diffusion, we deduce how diffusion affects the dynamics of the 
adsorbate configuration during the desorption process. 

3 Initial Configurations 

Figure [T] shows typical snapshots of the initial configurations that we used for the des- 
orption simulations. Each of these typical configurations yields a correlation length, 
obtained from averaging 100 statistically similar initial configurations. Each simulation 
setup is identified by this initial correlation length. Figure Q] illustrates the fact that 
within the constraint of our simulation setup, larger correlation length at a given cov- 
erage results in a cluster size distribution dominated by larger clusters. Figure [2] shows 
number density histograms for the initial configurations. These figures show that we 
have four initial configurations with size distributions that span (a) 1 < s < 12000 for 
f init = 4.96, (b) 1 < s < 3000 for £ Mt = 2.93, (c) 1 < s < 1500 for £ Mt = 2.15, and (d) 
1 < s < 400 for Cinit = 1-55. 

4 Simulation Results 

As the adsorbate configurations changed during the course of the desorption process, 
we measured the coverage without and with diffusion. Figure [3] shows 6(t) produced by 
desorption starting from each of the initial configurations. As expected, a larger initial 
correlation length results in slower desorption as the initial configuration is dominated 
by larger clusters. 

To quantify the effect of diffusion on 6(t), we define Ati/ 2 = ti/2 ~ ^'xit Here the 
half-time i]/ 2 is defined as the time for the system to reach a coverage of #mit/2 = 0.175 
without diffusion, and t^ 2 is the half-time with diffusion. We used Ati/ 2 to measure the 
acceleration effect due to the diffusion. Negative Ati/ 2 , therefore, indicates retardation 
of the desorption process. 

To measure the effect of diffusion on the cluster size distribution p(s), we compared 
the changes without and with diffusion at a given 9 by measuring and p(s). FigureH] 
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shows for each of the four initial configurations. To quantify the effect of diffusion 
on £(0), we define (S) = j^- Y^q {£,'(0) ~ £(^))> where £(0) and £'(0) are the correlation 
lengths at a given for desorption without and with diffusion, respectively, and is 
the number of terms in the sum. From Fig. Slit is evident that diffusion always increases 
the correlation length £ at a given 0. Figure [5] shows the results for Ai]/ 2 and (5) for 
1.5 < £ < 8.1. The inset in Fig. [5](a) shows the cross-over to retardation (Ai^ < 0) 
for small £i n it- 

To further explore the changes in the cluster size distribution caused by diffusion, 
we sampled p(s) histograms with and without diffusion at = 0.25, as shown in Fig. [6l 
For £ = 4.96 (fig. EJ^a)), there are no significant changes by diffusion, indicating that 
diffusion in this case only accelerates the desorption, without having any significant 
effect on the size distribution. In Fig. [6jb,c,d) we see depletion of small clusters and 
enhancement of large clusters, indicating coarsening due to diffusion. 

5 Discussion 

From the results for 0{t) (Fig. [3]) we find that when starting with larger £i n it> the 
acceleration effect of diffusion (as shown by Ai]/ 2 ) increases. This effect decreases 
when we start with smaller £i n it- At a certain point, the effect crosses over to slowing 
down the desorption (Ati/ 2 becomes negative for the smallest £j n i t (see FigJS^a)). The 
measurement of £(0) (Figs. I4)5h . on the other hand, shows that (5) is always positive, 
and the magnitude of the change increases as we start with smaller £i n it (see Fig. E^b)). 

Figure E^b,c,d) shows depletion of small clusters by diffusion. We associate this 
depletion and the positive (5) with a coarsening effect. The small positive (5) for 
the simulation with £i n it = 4.96 indicates a small coarsening effect. In this particular 
simulation, the acceleration effect dominates. For simulations with £i n j t = 2.93 and 
Cinit = 2.15, we observe larger (5) and smaller Ati/ 2 - Finally, for £i n j t = 1.55, we see 
large (5) with negative Ai]/ 2 , which means that the coarsening effect dominates in this 
case. 

In conclusion, lateral diffusion of adparticles creates two competing effects during 
desorption: an acceleration effect and a coarsening effect. Increasing £; n i t decreases the 
coarsening effect and increases the acceleration effect, due to enhanced release of par- 
ticles from the perimeters of large clusters. Further details will be presented elsewhere 
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£ = 4.96 £ = 2.93 £ = 2.15 £ = 1.55 




Figure 1: Typical initial snapshots prepared by adsorption at different chemical poten- 
tials, \i — fio = 0.4, 1.0, 2.0, and 10.0 from left to right. 
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Figure 2: Number density plots p(s) for initial configurations prepared by adsorption, 
similar to the ones shown in Fig. [TJ 
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Figure 3: The effects of increased diffusion on 0(t), starting with #; n it = 0.35 and 
different £i n i t . Note the different time scales. 
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Figure 4: (Color online) Correlation length £ as a function of coverage 9, measured 
without and with diffusion. This plot shows that diffusion always has a coarsening 
effect of varying magnitude, depending on the initial correlation length £; n it- Note that 
6 decreases from left to right in this figure, as the time increases. 
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Figure 5: Diffusion effects on 0(t) and £(0), obtained by measuring (a) Ati/ 2 and (b) 
(5). (See definitions in the text.) Inset in (a): crossover to negative Ai]/ 2 . 
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Figure 6: Size distributions for 9 = 0.25 without and with diffusion. The simulations 
start from 9- m it = 0.35, with (a) £i n it = 4.96, (b) £i n it = 2.93, (c) £i n i t = 2.15, and (d) 
6nit = 1-55. Note the different scales on both axes in the different parts of this figure. 
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